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User's  Guide  and  Documentation 
for  the  NRL  Environmental  Data 
Browsing  System  (NEBS) 

1.  Introduction.  The  NRL  Environmental  data  Browsing  System  (NEBS)  was  developed  to 
quickly  visualize  4D  (x,  y,  z,  t)  gridded  field  data  from  numerical  models  and  to  help  evaluate 
and  diagnose  model  behavior.  The  NEBS  allows  the  user  to  quickly  slice  and  animate,  data 
displays  in  x,  y,  z,  or  t.  NEBS  was  specifically  developed  to  support  the  Variability  of  Coastal 
Atmospheric  Refiractivity  (VOCAR)  experiment,  but  was  designed  to  be  very  general.  NEBS  can 
be  used  to  browse  almost  any  data  formatted  in  accordance  with  the  specifications  in  Appendix 
A.  The  script  for  the  NEBS  is  listed  in  Appendix  B. 

2.  General  Features.  The  NEBS  uses  the  Grid  Analysis  and  Display  System  (GrADS;  Doty  and 
Kinter,  1992)  developed  by  the  Center  for  Ocean-Land-Atmosphere  Interactions  (COLA)  at  the 
University  of  Maryland.  NEBS  is  written  in  the  GrADS  scripting  language  which  GrADS 
interprets  and  executes.  GrADS  is  supported  on  a  variety  of  platforms  including  DOS-based  PC 
and  Unix  workstations  (Digital,  IBM,  Hewlett-Packard,  NeXT,  Silicon  Graphics,  Sun,  and 
others).  However,  the  NEBS  is  currently  only  supported  on  the  Unix  workstations. 

2.1  Data  Graphical  User  Interface.  The  NEBS  is  a  graphical  interface  designed  for  quick-look 
and  diagnostics  of  model  output.  It  reads  data  in  the  GrADS  data  format  and  other  data  formats 
supported  by  GrADS,  such  as  GRIB.  Simple  translation  programs  are  available  to  reformat 
NORAPS  and  NOGAPS  files. 

2.2  Mouse  and  Button  Control.  Browsing  is  controlled  entirely  with  a  mouse  cursor,  mouse 
buttons,  and  GUI  "buttons"  to  select  options  and  execute  commands. 

2.3  Selectable  Variables.  A  button  for  each  variable  in  a  data  set  is  created  automatically  when  a 
data  file  is  selected.  Multiple  variables  can  be  selected  and  displayed  simultaneously  as  contours. 
One  variable,  at  most,  can  be  shaded. 

2.4  Selectable  Displays.  Vector  variables  can  be  displayed  as  streamlines,  vectors,  or  barbs. 
Contour  intervals  and  limits  can  be  customized  for  each  variable,  and  vector  variables  can  be 
parsed,  if  too  dense. 

2.5  Selectable  Dimensions  and  Times.  Spatial  dimension  and  time  bars  appropriate  for  the  data 
contents  of  the  selected  file  are  displayed  and  allow  the  selection  of  precise  levels  or  varying 
ranges  in  each  dimension  or  time  (only  three  of  these  four  parameters  can  vary  simultaneously). 

2.6  Selectable  Stills  or  Animation.  Maps,  cross  sections,  or  line  plots  are  displayed,  or  animated, 
based  on  the  selected  dimensional  conditions.  This  allows  the  data  to  be  sliced  and  animated  in 
X,  Y,  Z,  or  T.  Animation  occurs  if  the  indicated  dimension  or  time  bars  have  highlighted  ranges. 

2.7  Time  Interpolation.  Time  animations  can  be  more  highly  resolved  by  selecting  a  time 
interval  to  interpolate  the  data  set  values. 

2.8  Zoom  and  Pan.  Zoom  in,  zoom  out,  and  panning  of  map  displays  of  model  data  can  be  used 
for  a  quick  investigation  of  data  anomalies. 
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2.9  Selectable  Reference  Observations.  Radiosonde  profiles  and  surface  observations  can  be 
compared  to  model  forecasts  for  evaluation. 

Profiles.  Raob  data  and  model  profiles  can  be  displayed  simultaneously  for  only  a  single 
station.  The  profiles  can  be  time  animated,  but  if  time  interpolated,  the  raob  data  only  appears  at 
its  valid  time. 

Time  Lines.  Maps  of  selected  model  variables  and  observed  surface  winds  can  be 
displayed  and  animated.  If  time  interpolated  in  a  map  view,  both  model  variables  and  surface 
observations  are  interpolated  and  displayed  simultaneously.  Individual  surface  stations  can  be 
independently  selected  in  the  map  view  to  allow  the  display  of  selected  model  variables  and 
surface  observations  as  time  lines  for  each  station. 

2.10  Selectable  Output  Media.  Once  a  single  display  or  animation  sequence  is  specified,  the 
frame(s)  can  be  sent  to  a  black/white  printer,  color  printer,  GrADS  movie  file,  or  video  editor 
(not  currently  implemented)  for  creation  of  a  VHS  video  tape. 


3.  InstaHation  and  Setap. 

If  GrADS  is  not  already  installed,  it  will  be  necessary  to  download  Version  1.5,  or  later,  from 
COLA  at  the  University  of  Maryland  (grads.iges.com,  IP#  192.239.84.50). 

If  operating  remotely  over  a  LAN,  setup  the  local  xhost  to  accept  a  remote  display  using  the 
‘xhost  +’  command  on  your  local  workstation  and  designate  your  workstation  to  the  xserver 
using; 

setenv  DISPLAY  workstation:0.0 

The  directory  on  the  remote  xserver  where  GrADS  is  installed  must  be  added  to  your  search 
PATH  and  the  environment  must  be  set  as  follows  in  your  login  file: 

setenv  TERM  xhost 
setenv  GADDIR  directory 

The  GrADS  data  control  files  for  the  data  to  be  browsed  must  be  in  the  same  directory  as  the 
NEBS  script.  If  station  raob  and  surface  data  files  are  used,  they  must  be  accompanied  by  one 
file  each  that  specifies  the  station  names  and  locations  in  each  observed  data  file;  named 
’stations.ua'  and  'stations.sfc',  respectively.  The  formats  for  these  files  are  described  in  Appendix 
A.  This  naming  convention  is  necessary  unless  the  NEBS  script  is  modified.  GrADS  control 
files  for  model  data  may  have  any  name,  but  must  have  the  standard  “.ctl”  extension. 


4.  Using  the  NRL  Environmental  Browslmg  System. 

4.1  Start  Up.  First,  open  an  xterm  window,  change  directories  to  the  directory  containing  the 
browser,  and  initiate  GrADS  by  typing 
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grads  -1  {Enter}. 


At  the  GrADS  prompt  (i.e.,  ga>),  type 

run  nebs.gs  {Enter}. 

Alternatively,  these  can  be  combined  into  a  single  command  in  the  xterm  window  by  typing 

grads  -Ic  “run  nebs.gs”  {Enter}. 

A  list  of  available  data  files  will  appear.  A  file  from  this  list  may  be  typed  in  response  to  the  query 
for  a  grid  data  filename.  (This  may  also  be  accomplished  with  'copy'  and  'paste'  operations,  if 
these  features  are  available  on  your  workstation.)  Next,  a  query  will  appear  requesting  the  surface 
data  filename,  followed  by  a  query  for  the  raob  data  filename.  If  surface  or  raob  data  are  not 
available,  or  are  not  desired,  an  {Enter}  will  bypass  these  queries.  The  browser  will  start  up,  as 
shown  in  Figure  1,  by  displaying  a  window  with  control  buttons  for  selecting  variables,  graph 
types,  time  settings,  dimension  settings,  and  basic  functions,  such  as  reset,  print,  quit,  and  display. 


Figure  1.  Basic  NEBS  control  buttons. 
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4.2  Variable  Selection.  Variables  can  be  "toggled"  off  and  on  by  placing  the  mouse  cursor  on  the 
desired  variable's  button  and  by  clicking  the  left  or  middle  mouse  button.  Also,  the  browser  will 
start  automatically  with  sea  level  pressure  or  surface  pressure  displayed  (Fig.  2),  if  present  in  the 
data  set,  to  provide  an  initial  orientation  to  the  meteorological  conditions  existing  at  the  first  time 
step. 


Figure  2.  Initial  pressure  display  (when  data  present). 

4.3  Variable  Customization.  A  menu  of  options  (Fig.  3)  can  be  obtained  by  placing  the  mouse 
cursor  on  the  desired  variable's  button  and  by  clicking  the  right  mouse  button.  All  options  must  be 
selected  using  the  right  mouse  button.  Clicking  the  left  or  middle  mouse  button  accepts  the 
displayed  menu  values. 

4.3.1  Contour  labels.  The  presence  of  contour  labels  can  be  toggled  on  and  off  by  clicking 
the  On/Off  status  button. 


4.3.2  Contour  interval.  The  plot  contour  (or  shading) 
decremented  by  clicking  the  ’>'  or  '<’  symbol,  respectively. 


interval  can  be  incremented  or 
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Labels 

pmsi  interval 
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LINE 
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Marker  none 


Figure  3.  Variable  customization  menu. 

4.3.3  Contour  maximum.  The  maximum  value  contoured  (or  shaded)  can  be  incremented 
or  decremented  by  clicking  the  '>'  or  '<’  symbol,  respectively. 

4.3.4  Contour  minimum.  The  minimum  value  contoured  (or  shaded)  can  be  incremented 
or  decremented  by  clicking  the  ’>'  or  '<'  symbol,  respectively. 

4.3.5  Line  style.  The  style  of  line  used  for  plotting  can  be  selected  as  ‘solid’,  ‘long  dash’, 
‘short  dash’,  or  ‘dotted’  by  clicking  on  the  ‘?’  symbol  until  the  desired  option  appears. 

4.3.6  Line  marker.  The  marker  for  plotting  data  can  be  selected  as  ‘none’,  ‘cross’,  ‘open 
circle’,  ‘closed  circle’,  ‘open  square’,  or  ‘closed  square’  by  clicking  on  the  ‘?’  symbol  until  the 
desired  option  appears. 

4.3.7  Skip  interval.  Vector  variables  can  be  decimated  by  skipping  grid  points.  When  skip 
is  set  to  'O'  every  point  is  plotted.  If  skip  is  set  to  T,  '2',  or  'n',  the  2nd,  3rd,  or  n+lth  points  in 
latitude  and  longitude  are  plotted. 
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4.4  Graph  Type  Selection.  By  default,  the  graph  type  that  is  highlighted  when  a  variable  is 
selected  will  be  used  to  plot  that  variable. 

4.4. 1  Contour.  This  graph  type  is  the  most  common  for  scalar  variables  and  can  be  for  as 
many  variables  as  desired.  Each  subsequent  usage  of  contouring  within  the  same  graph  will 
produce  a  different  contour  color  and  corresponding  label  color. 

4.4.2  Shaded.  The  'shaded'  graph  type  can  only  be  used  once  for  a  scalar  variable  since  a 
second  selection  would  overdraw  the  first  plotted  variable.  If  selected  a  second  time,  a  warning  is 
issued  and  the  name  of  the  first  shaded  variable  is  given  so  its  graph  type  may  be  changed  or  the 
variable  may  be  deselected. 

4.4.3  Vector.  This  graph  type  places  vectors  with  their  tail  at  their  origin  plotted  in  the 
direction  the  wind  is  going.  It  is  primarily  for  wind  data  and  is  similar  to  the  shaded  graph  type  in 
that  it  can  only  be  used  for  one  pair  of  vector  variables. 

4.4.4  Barb.  This  graph  t5q)e  places  barbs  with  their  head  at  their  origin  plotted  in  the 
direction  from  where  the  wind  is  coming.  It  is  primarily  for  wind  data  and  is  similar  to  the  shaded 
graph  type  in  that  it  can  only  be  used  for  one  pair  of  vector  variables. 

4.4.5  Stream.  This  graph  type  produces  magnitude  colored  streamhnes  and  is  primarily  for 
wind  data  and  is  similar  to  the  shaded  graph  type  in  that  it  can  only  be  used  for  one  pair  of  vector 
variables. 

4.4.6  Mixed.  Contours  of  isotachs  can  be  overlaid  on  vector,  barb,  or  streamline  displays 
of  winds  by  selecting  one  wind  component  (i.e.,  u  or  v)  with  'contour'  and  selecting  the  other  wind 
component  with  'vector',  'barb'  or  'stream'. 

4.5  Date-Time  Selection.  The  time  span  contained  within  the  data  set  is  displayed  as  a  vertical  bar 
with  a  tick  mark  for  each  time  group  in  the  data.  Either  data  set  time  marks  or  interpolated  time 
marks  can  be  selected. 

4.5.1  Single  time.  A  single  time  value  is  selected  by  clicking  with  the  left  mouse  button  on 
the  desired  time  tick.  The  exact  time  and  date  of  the  selected  tick  is  displayed. 

4.5.2  Time  range.  The  first  time  value  is  selected  using  either  the  left  or  middle  mouse 
button.  Then  the  second  time  value  is  selected  with  the  middle  mouse  button.  The  selected  range 
is  highlighted  and  the  beginning  and  ending  date-time  values  are  displayed. 

4.5.3  Time  settings.  The  current  time  interval  between  tick  marks  is  displayed  vertically  at 
the  center  of  the  time  bar.  Clicking  with  the  right  mouse  button  on  the  time  interval  displays  the 
time  settings  menu  (Fig.  4).  Options  for  displaying  a  local  clock  or  setting  the  desired  time 
interpolation  interval  are  available.  Options  must  be  selected  using  the  right  mouse  button. 
Clicking  the  left  or  middle  mouse  button  accepts  the  displayed  time  option  settings. 


6 


Figure  4.  Time  settings  menu. 


4.5.3. 1  Local  time  clock.  A  clock  face  showing  the  local  time  can  be  toggled  on  or 
off  by  clicking  its  On/Off  status  button.  The  local  time  offset  is  the  hour  calculated  using  the 
display’s  center  longitude  rounded  to  the  nearest  15°  increment  from  GMT. 

4.5. 3.2  Time  interpolation.  The  time  interval  can  be  increased  or  decreased  by 
clicking  the  ’>'  or  '<'  symbol,  respectively.  If  an  interpolated  interval  is  selected  the  vertical  time 
bar  indicates  the  status  (Fig.  5)  and  the  interpolated  time  intervals  are  displayed  with  yellow  tick 
markers  on  the  time  bar. 

4.6  Dimension  Selection.  The  default  spatial  dimension  is  level  (Z)  and  may  be  quoted  in  feet, 
meters,  sigma  levels,  etc.,  as  determined  by  the  data  set.  Longitude  (X)  or  latitude  (Y)  in  degrees 
may  be  selected  by  clicking  on  the  appropriate  button  using  any  mouse  button.  The  range  and 
selected  dimensional  values  are  displayed  on  a  vertical  bar  similar  to  time. 

4.7  Level  Selection.  Select  'Lev'  using  any  mouse  button.  The  vertical  extent  of  the  data  set  is 
displayed  as  a  vertical  bar  with  a  tick  mark  for  each  available  vertical  level. 
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Figure  5.  Customization  of  time  interval. 


4.7.1  Single  level.  A  single  level  value  is  selected  by  clicking  with  the  left  mouse  button 
on  the  desired  level's  tick  mark.  The  exact  level  of  the  selected  tick  is  displayed. 


4.7.2  Level  range.  The  first  level  value  is  selected  using  either  the  left  or  middle  mouse 
button.  Then  the  second  level  value  is  selected  using  the  middle  mouse  button.  The  selected 
range  is  highlighted  and  the  beginning  and  ending  level  values  are  displayed. 

4.8  Latitude  Selection.  Select  'Lat'  using  any  mouse  button.  The  north-south  extent  of  the  data  set 
is  displayed  as  a  vertical  bar  with  a  tick  mark  for  each  available  latitude. 

4.8.1  Single  latitude.  A  single  latitude  value  is  selected  by  chcking  with  the  left  mouse 
button  on  the  desired  latitude's  tick  mark.  The  exact  latitude  of  the  selected  tick  is  displayed. 


4.8.2  Latitude  range.  The  first  latitude  value  is  selected  using  either  the  left  or  middle 
mouse  button.  Then  the  second  latitude  value  is  selected  with  the  middle  mouse  button.  The 
selected  range  is  highlighted  and  the  beginning  and  ending  latitude  values  are  displayed. 
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4.9  Longitude  Selection.  Select  'Lon'  using  any  mouse  button.  The  east-west  extent  of  the  data  set 
is  displayed  as  a  vertical  bar  with  a  tick  mark  for  each  available  longitude. 

4.9. 1  Single  longitude.  A  single  longitude  value  is  selected  by  clicking  with  the  left  mouse 
button  on  the  desired  longitude's  tick  mark.  The  exact  longitude  of  the  selected  tick  is  displayed. 

4.9.2  Longitude  range.  The  first  longitude  value  is  selected  using  either  the  left  or  middle 
mouse  button.  Then  the  second  longitude  value  is  selected  with  the  middle  mouse  button.  The 
selected  range  is  highlighted  and  the  beginning  and  ending  longitude  values  are  displayed. 

4.10  Display  and  Animation.  The  selected  graph  types,  variables  and  conditions  are  displayed  by 
clicking  on  'DISPLAY'  with  any  mouse  button,  or  by  clicking  the  right  mouse  button  in  the  central 
clear  area  of  the  window.  (Note:  Clicking  on  a  dimension  bar  or  time  bar  with  the  right  mouse 
button  win  select  a  new  single  dimension  or  time  value  and  immediately  re-display  the  selected 
variables  and  graph  types.)  Display  and  animation  actions  are  inferred  from  the  fixed  level,  or 
varying  range  status  of  the  time  bar  and  the  indicated  dimension  bar. 

4.10.1  Single  frame.  In  general,  if  a  single  fixed  time  is  selected  and  the  indicated 
dimension  has  a  single  fixed  value  selected,  only  a  single  frame  will  be  displayed.  Exceptions  are 
shown  in  Table  1  and  wiU  be  discussed  in  Section  5. 

4.10.2  Animation.  In  general,  if  the  time  bar  shows  a  varying  range  selected,  an  animation 
in  time  will  result.  In  general,  if  the  indicated  dimension  bar  shows  a  varying  range  is  selected,  an 
animation  in  that  dimension  will  result.  Exceptions  are  shown  in  Table  1  and  will  be  discussed  in 
Section  5. 


Table  1 .  Gridded  data  display  options. 


5.  Plotting  Gridded  Data. 

5.1  Plot  Types.  Various  combinations  of  fixed  and  varying  conditions  for  latitude,  longitude, 
level,  and  time  can  be  selected.  These  combinations  further  interact  with  whether  the  indicated 
dimension  is  fixed  or  varying.  The  resulting  plot  types  for  the  combinations  shown  in  Table  1 
are  discussed  below: 

5.2  No  Plot. 

-  All  dimensions  and  time  are  fixed.  The  values  of  the  selected  dimensions,  time,  and 
variables  are  printed  (Fig.  6). 

-  All  dimensions  and  time  are  varying.  This  case  is  prohibited. 


Figure  6.  Printed  data  for  fixed  conditions. 
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5.3  Line  Plot. 

-  One  dimension  or  time  is  varying.  Line  plots  of  the  selected  variables  are  displayed  in 
different  colors  with  color  coded  axes  and  labels  (Fig.  7). 

-  Time  and  a  hidden  dimension  are  varying.  A  time  animation  of  a  line  plot  for  the 
hidden  dimension  will  occur. 

-  The  indicated  dimension  and  a  hidden  dimension  are  varying.  A  slice  animation  of  a 
line  plot  for  the  hidden  dimension  will  occur  with  the  indicated  dimension  varying. 


nor^ps  fields  :i 
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Figure  7.  Typical  line  plot. 

5.4  2D  Section. 

-  One  dimension  and  time  are  varying  (Fig.  8).  A  static  time  cross  section  will  occur  for 
the  indicated  dimension. 

-  Two  dimensions  are  varying,  but  time  is  fixed  (Fig.  2).  Map  plots  with  coast  lines  and 
political  boundaries  are  displayed  if  level  is  fixed.  Vertical  cross  sections  are  displayed  along  the 
selected  latitude  (longitude),  if  latitude  (longitude)  is  fixed  (Fig.  9). 
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Figure  8.  Typical  Z-T  time  section. 

-  Two  dimensions  are  varying  and  time  is  varying.  A  time  animation  of  a  map  or  cross 
section  will  occur  if  the  indicated  dimension  is  fixed,  or  time  section  will  animate  in  the 
indicated  dimension  if  it  is  varying. 

-  Three  dimensions  are  varying.  A  slice  animation  of  the  maps  will  occur  if  level  is  the 
indicated  dimension,  or  a  slice  animation  of  vertical  cross  sections  will  occur  for  latitude 
(longitude)  if  latitude  (longitude)  is  the  indicated  dimension. 

6.  Plotting  Surface  Data. 

6. 1  Surface  Station  Availability.  If  surface  station  data  are  available  for  comparison  with  model 
surface  forecasts,  a  'SURFACE'  button  will  appear  under  the  row  of  variable  buttons.  Clicking 
on  'SURFACE'  with  any  mouse  button  will  toggle  the  'SURFACE'  button  on  and  off. 
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Figure  9.  Typical  Z-X  vertical  cross  section 


6.2  Surface  Station  Wind  Display.  Surface  station  wind  data  are  displayed  when  a  map  display  is 
selected  (i.e.,  latitude  and  longitude  varying  with  level  and/or  time  fixed).  If  a  time  range  is 
selected,  the  selected  variable  and  the  surface  wind  barbs  will  animate.  If  a  time  interpolation 
interval  is  selected,  both  the  model  data  and  the  surface  wind  data  will  be  interpolated  and 
displayed.  If  a  range  of  levels  is  selected,  the  model  data  will  animate  in  Z,  but  the  surface  winds 
will  remain  fixed. 

6.3  Surface  Station  Selection.  Surface  stations  can  be  selected  only  from  a  fixed  X-Y  map 
display,  without  any  animation  in  Z  or  T  because  the  selection  and  animation  processes  are 
incompatible.  The  X-Y  map  display  (Fig.  10)  of  any  variable  will  be  overlaid  with  surface  wind 
barbs  for  all  available  surface  stations.  When  station  markers  are  visible,  individual  stations  can 
be  selected  by  clicking  on  the  desired  station's  marker  using  the  middle  mouse  button.  The 
station  will  be  highlighted  and  its  3 -character  nomenclature  will  be  displayed  to  its  upper  right. 
Repeated  clicking  on  a  station  with  the  middle  mouse  will  toggle  its  selection  on  and  off. 
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Figure  10.  Display  for  surface  station  selection. 

6.4  Surface  Time  Line  Display.  When  the  desired  surface  stations  have  been  selected,  clicking 
the  middle  mouse  button  outside  the  X-Y  map  display  will  cause  the  display  of  stacked  time 
lines,  one  for  each  selected  surface  station  (Fig.  1 1).  Variable  selections  can  be  changed  and  new 
station  time  lines  can  be  displayed  by  clicking  on  new  variables  followed  by  an  immediate  click 
with  the  middle  mouse  button  in  the  central  clear  area  of  the  window. 


15 


Figure  11.  Typical  surface  time  lines. 


6.5  Model  Time  Section  Display.  Once  stations  have  been  selected,  it  is  possible  to  plot  a  Z-T 
section  of  model  data  at  the  selected  location  by  first  defining  the  desired  range  in  Z  by  using  the 
middle  mouse  button  on  the  level  bar  and  then  clicking  in  the  center  of  the  window  with  the 
middle  mouse  button  (Fig.  12).  (Note:  To  change  the  station  selection  you  must  first  reset  the 
level  to  a  single  value  of  Z  or  the  display  will  animate  in  Z  when  the  map  display  is  activated.) 

6.6  Surface  Display  Customization.  The  contour  limit  selections  available  under  the  variable 
menus  can  be  changed  to  control  the  plot  range  of  the  selected  variables.  If  no  time  range  is 
selected  the  default  is  the  full  time  range.  Selecting  a  time  range  will  expand  the  display  in  time 
to  the  desired  interval.  Model  time  section  variables  can  be  shaded  or  displayed  as  vectors. 
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Figure  12.  Typical  model  time  section. 

7.  Plotting  Upper  Air  Data. 

7.1  Upper  Air  Station  Availability.  If  upper  air  station  data  are  available  for  comparison  with 
model  profiles,  a  row  of  station  buttons  will  appear  under  the  row  of  variable  buttons  with  the  3- 
character  nomenclature  of  each  station  indicated. 

7.2  Upper  Air  Station  Selection.  Clicking  on  a  station's  button  with  any  mouse  button  will  toggle 
that  station  on  and  off.  Since  upper  air  data  can  be  displayed  for  only  one  station,  any  previously 
selected  stations  will  be  deselected. 

7.3  Upper  Air  Data  Selection.  Variables  aro  selected  the  same  as  gridded  and  surface  data  by 
using  the  left  or  middle  mouse  button  to  toggle  the  variable  selection. 

7.4  Upper  Air  Station  Data  Display.  Upper  air  data  and  model  vertical  profiles  are  displayed 
(Fig.  13)  the  same  as  gridded  data  by  either  clicking  ’DISPLAY"  or  clicking  the  right  mouse 
button  in  the  window's  central  area. 
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Figure  13.  Typical  radiosonde  and  model  profiles. 


8.  Zooming  and  Panning  Map  Displays. 

Only  X-Y  map  displays  of  a  gridded  data  variable  can  be  zoomed  in,  zoomed  out,  and  panned. 

8.1  Zooming  In.  The  first  click  within  a  map  display  using  the  left  mouse  button  will  cause  the 
scale  to  be  halved  and  the  display  to  be  re-centered  on  the  location  of  the  mouse  cursor.  A 
subsequent  left  mouse  button  click  will  zoom  in  only  if  near  the  map  center.  Each  zoom  results 
in  a  halving  of  the  map  scale,  or  a  size  doubling  of  map  features. 

8.2  Panning.  If  the  left  mouse  button  click  is  outside  the  central  20%  of  the  map,  the  scale  will 
remain  the  same,  but  the  map  will  be  re-centered  on  the  new  mouse  cursor  location.  This  allows 
the  ability  to  pan  without  further  zooming  until  the  feature  of  interest  is  near  the  center,  when 
zooming  in  can  be  continued  if  desired. 
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8.3  Zooming  Out.  The  assigned  plot  area  has  a  hidden  boundary  that  is  sensitive  to  clicks  with 
the  left  mouse  button.  Generally,  a  left  mouse  button  click  on  or  outside  the  map  border  will 
result  in  a  doubling  of  the  map  scale.  The  area  under  the  map  has  a  larger  sensitive  band  that 
provides  a  more  reliable  zoom  out  activation. 


9.  Printing  Displays  and  Animation  Sequences. 

Executing  any  print  option  re-displays  the  last  display  request  and  sends  it  to  the  selected 
destination.  Clicking  with  any  mouse  button  on  'PRINT'  will  display  the  print  menu  (Fig.  14). 
Clicking  on  a  selection  will  activate  that  choice.  A  large  number  of  frames  can  be  easily 
requested  which  may  overload  the  selected  printer  and  local  disk  storage  capacity.  The  number 
of  frames  is  quoted  and  a  final  choice  to  continue  is  required.  The  print  function  is  implemented 
in  a  companion  script  called  ‘output.gs’  which  must  be  edited  to  include  print  commands 
appropriate  to  your  system. 


Figure  14.  Print  menu. 
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9.1  BW.PLOT.  The  requested  frames  are  plotted  to  a  local  black  and  white  printer,  or  a  network 
printer,  as  defined  in  the  ‘output.gs’  script. 

9.2  C.PLOT.  The  requested  frames  are  plotted  to  a  network  color  printer,  as  defined  in  the 
‘output.gs’  script. 

9.3  MOVIE.  The  requested  frames  are  plotted  in  a  frame  animator  on  the  host  display.  This 
allows  a  quick  preview  of  the  requested  frames  before  they  are  printed  or  sent  to  the  video  editor. 

9.4  VIDEO.  The  requested  frames  are  properly  reformatted  and  sent  to  a  commercial  video 
editor  where  title  frames  can  be  added  and  the  frames  can  be  recorded  to  VHS  video  tape.  (This 
currently  is  not  implemented.  The  necessary  commands  and  data  paths  must  be  set  up  in  the 
companion  GrADS  script  ‘outputgs’.) 

10.  Resetting. 

Occasionally,  it  may  be  convenient  to  return  to  the  initial  settings  that  existed  when  NEBS  was 
started.  Rather  than  resetting  many  conditions  individually,  any  mouse  button  click  on  the 
'RESET'  button  will  re-initiahze  all  variables  and  conditions. 

11.  Exiting. 

The  NEBS  session  should  be  terminated  by  a  mouse  button  click  on  the  'QUIT'  button.  This 
returns  control  to  the  GrADS  prompt  in  the  originating  xterm  window. 


The  NRL  Environmental  data  Browsing  System  (NEBS)  is  able  to  browse  any  data  set  provided 
that  it  is  in  a  GrADS  compatible  data  set  format.  NEBS  is  capable  of  easily  producing 
animations  in  time  or  a  selected  dimension,  x,  y,  or  z.  These  animations  allow  a  quahtative 
assessment  of  data  consistency  and  model  behavior  and  easily  reveal  anomalies  and  erroneous 
trends.  By  setting  time  and  spatial  dimensions  at  the  location  of  an  identified  problem,  NEBS 
gives  a  quantitative  listing  of  individual  variables.  Variables  and  dimensions  can  be  easily  set 
using  mouse  button  clicks  on  the  displayed  GUI  buttons.  Graphics  are  displayed  quickly  and  can 
be  animated  on  most  host  platforms  at  about  one  frame  per  second.  Over  a  high  speed  LAN  this 
degrades  to  about  one  frame  per  two  seconds.  Selected  data  displays  can  be  easily  printed,  in 
black  &  white  or  color,  or  converted  to  on-screen  animations  or  video  VHS  tape. 


GrADS  provides  excellent  data  handling  and  2-D  graphical  display  of  gridded  data.  However, 
upgrading  GrADS  or  coupling  GrADS  to  a  3-D  graphics  system,  such  as  VIS-5D,  is  highly 
desirable  to  provide  3-D  volume  renderings  of  isosurfaces  and  wind  fields.  NEBS  should  be 
revised  whenever  a  3-D  rendering  capability  becomes  available  for  GrADS. 


20 


ACBCNOWLEDGEMENTS 


The  original  concept  for  the  Naval  Environmental  data  Browsing  System  was  inspired  by  the 
GrADS  data  browser  created  by  Steve  J.  Lord  at  the  National  Meteorological  Center.  The 
GrADS  author,  Brian  Doty,  and  other  GrADS  users,  especially.  Dr.  Michael  Fiorino  were  very 
helpful  in  the  NEBS  development.  Also  Dr.  Stephen  Burk  and  Dr.  William  Thompson  of  NRL 
are  thanked  for  their  guidance  and  feedback  during  NEBS  development.  The  support  of  the 
sponsor,  the  Space  and  Naval  Warfare  Systems  Command  under  program  element  0603207N,  is 
gratefully  acknowledged. 


REFERENCES 


Doty,  B.  and  J.  Kinter,  1992:  The  Grid  Analysis  and  Display  System  (GrADS):  A  practical  tool 
for  earth  science  visualization.  Eighth  International  Conference  on  Interactive  Processing 
Systems  for  Meteorology,  Oceanography  and  Hydrology,  Jan  1992,  Atlanta,  GA  (Amer.  Meteor. 
Soc.),  115-116. 


Doty,  B.  and  J.  Kinter,  1993:  The  Grid  Analysis  and  Display  System  (GrADS):  An  Update, 
Ninth  International  Conference  on  Interactive  Processing  Systems  for  Meteorology, 
Oceanography  and  Hydrology,  Jan  1993,  Anaheim,  CA  (Amer.  Meteor.  Soc.),  165-167. 


21 


This  Page  Intentionally  Left  Blank. 


22 


Appendix  A:  Data  Input  Formats 


A.l.  Gridded  Data  Sets . A-3 

A.2.  Station  Data  Sets . A-12 

A.3.  Examples  of  Creating  a  Gridded  Data  Set . A- 16 

A.4.  Examples  of  Creating  Station  Data  Sets . A- 17 

A.5.  Examples  of  ‘stations.ua’  and  ‘stations.sfc’  Data  Sets  .  A-20 


A-l 


This  Page  Intentionally  Left  Blank. 


A-2 


GrADS  Data  Sets’ 

A.  1 .  Gridded  Data  Sets. 

The  GrADS  gridded  data  set  is  a  direct  access  binary  data  set.  It  may  contain  any  number  of 
variables  at  specified  longitude,  latitude,  vertical,  and  time  intervals. 

GrADS  views  this  data  set  as  a  giant  array  -  with  X  (longitude)  varying  the  fastest,  then  Y 
(latitude),  then  Z  (vertical  level),  then  the  variable  type,  then  T  (time). 

It  is  easier  for  us  to  think  of  the  data  set  in  terms  of  a  sequence  of  horizontal  grids,  where 
longitude  and  latitude  vary.  Each  horizontal  grid  represents  a  particular  variable  at  a  particular 
height  and  time.  Each  horizontal  grid  is  the  same  size  in  any  particular  GrADS  data  set  (if  you 
have  grids  of  different  sizes,  you  must  create  separate  data  sets). 

These  grids  are  written  to  the  data  set  in  the  following  order:  starting  with  a  particular  variable, 
grids  for  each  vertical  level  (at  a  particular  time)  are  written  out  in  ascending  order.  Then  the 
grids  for  the  next  variable  are  written  out.  When  all  the  grids  at  a  particular  time  have  been 
written,  grids  for  the  next  time  are  written. 

The  format  of  this  data  set  is  thus  exactly  the  same  as  the  COLA  Pressure  History  format,  except: 
there  are  no  date/time  records,  and  latitude  varies  from  south  to  north  (not  north  to  south  as  in  the 
pressure  history  data). 

Each  binary  gridded  data  set  is  described  by  a  data  descriptor  file,  essentially  a  table  of  contents 
for  the  binary  data  set.  Following  is  an  example  of  such  a  file: 

DSET  ua.dat 

TITLE  Upper  Air  Data 

UNDEF  -9.99E33 

XDEF  80  LINEAR  -140.0  1.0 

YDEF  50  LINEAR  20.0  1.0 

ZDEF  10  LEVELS  1000  850  700  500  400  300  250  200  150  100 
TDEF  4  LINEAR  0Z10aprl991  12hr 
VARS  5 

sip  0  0  sea  level  pressure 

z  10  0  heights 

t  10  0  temps 

td  6  0  dewpoints 

u  10  0  u  winds 

V  10  0  V  winds 

ENDVARS 


’  Extracted  from  GrADS  documentation  with  permission  from  Brian  Doty,  Center  for 
Ocean-Land- Atmosphere  Interactions  (COLA),  Department  of  Meteorology,  University  of  Maryland, 


The  data  descriptor  file  is  'free  foimat',  i.e.  each  entry  is  blank  dehmited  and  may  appear  in  any 
column.  Comment  records  start  with  an  asterisk  ('*')  in  column  1.  Comments  may  not  appear  in 
the  list  of  variable  records  (between  the  vars  and  endvars  records).  Records  may  not  be  more 
than  80  characters  long. 

In  this  example,  the  binary  data  set  is  named  ua.dat,  the  undefined,  or  missing,  data  value  is 
-9.99e33,  there  are  80  grid  points  in  the  X  direction,  50  in  the  Y  direction,  10  levels,  4  times,  and 
5  variables.  The  variables  z,  t,  u,  and  v  have  10  levels,  the  variable  td  has  6  levels,  and  the 
variable  sip  has  one  level  (see  below  for  a  more  specific  description  of  each  entry). 

Think  in  terms  of  the  X  and  Y  data  points  at  one  level  for  one  variable  at  one  time  being  a 
horizontal  grid.  This  grid  is  exactly  in  the  same  storage  order  as  a  FORTRAN  array,  in  this  case 
an  array  DIMENSION  A(80,50).  The  first  dimension  always  varies  from  west  to  east,  the  second 
from  south  to  north. 

In  the  above  example  the  horizontal  grids  would  be  written  in  the  following  order: 

Time  1 ,  Level  ?,  Variable  sip 
Time  1,  Level  1000,  Variable  z 
Time  1,  Level  850,  Variable  z 

then  levels  700,  500,  400,  300,  250,  200,  then 
Time  1,  Level  150,  Variable  z 
Time  1,  Level  100,  Variable  z 
Time  1,  Level  1000,  Variable  t 
Time  1,  Level  850,  Variable  t 

then  levels  700,  500,  400,  300,  250,  200,  then 
Time  1,  Level  150,  Variable  t 
Time  1,  Level  100,  Variable  t 
Time  1,  Level  1000,  Variable  td 
Time  1,  Level  850,  Variable  td 
Time  1,  Level  700,  Variable  td 
Time  1,  Level  500,  Variable  td 
Time  1,  Level  400,  Variable  td 
Time  1,  Level  300,  Variable  td 
Time  1,  Level  1000,  Variable  u 
Time  1,  Level  850,  Variable  u 

then  levels  700,  500, 400, 300, 250,  200,  then 
Time  i.  Level  150,  Variable u 
Time  1,  Level  100,  Variable  u 
Time  1,  Level  1000,  Variable  v 
Time  1,  Level  850,  Variable  v 

then  levels  700,  500, 400,  300, 250,  200,  then 
Time  1,  Level  150,  Variable  v 
Time  1,  Level  100,  Variable  v 
Time  2,  Level  ?,  Variable  sip 
Time  2,  Level  1000,  Variable  z 


A-4 


Time  2,  Level 
Time  2,  Level 
Time  2,  Level 
Time  2,  Level 


850,  Variable  z 
700,  Variable  z 
500,  Variable  z 
400,  Variable  z 


etc 

A  description  of  each  record  in  the  GrADS  data  descriptor  file  follows: 

DSET  data-set-name 

This  entry  specifies  the  name  of  the  binary  data  set.  It  may  be  entered  in  mixed 

case. 

If  the  binary  data  set  is  in  the  same  directory  as  the  data  descriptor  file,  you  may 
enter  the  filename  in  the  data  descriptor  file  without  a  full  path  name  by  prefixing  it 
with  a  ^  character.  For  example,  if  the  data  descriptor  file  is: 

/data/wx/grads/sa.ctl 

and  the  binary  data  file  is: 

/  data/wx/grads/sa.  dat 

you  could  use  the  following  file  name  in  the  data  descriptor  file: 

DSET  ^sa.dat 
instead  of: 

DSET  /data/wx/grads/sa.dat 

As  long  as  you  keep  the  two  files  together,  you  may  move  them  to  any  directory 
without  changing  the  entries  in  the  data  descriptor  file. 

TITLE  string 

A  brief  description  of  the  contents  of  the  data  set.  This  will  be  displayed  during  a 
QUERY  command,  so  it  is  helpful  to  put  meaningful  information  here. 

UNDEF  value 

The  undefined,  or  missing,  data  value.  GrADS  operations  and  graphics  routines 
will  ignore  data  with  this  value  from  this  data  set. 


A-5 


FILEHEADER  bytenum 

bytenum  -  the  number  of  bytes  of  non-data  at  the  beginning 
of  the  file 

Indicates  the  binary  data  file  has  a  header  where  bytenum  is  the 
number  of  bytes  in  the  header.  GrADS  will  skip  past  this  header, 
then  treat  the  file  as  though  it  were  a  normal  GrADS  file  after 
that  point.  This  option  is  valid  only  for  GrADS  gridded  data  sets. 

OPTIONS  <keywords> 

Some  keywords: 

options  <yrev>  <zrev>  <sequential>  <template> 

<big_endian>  <little_endian> 

yrev  -  y  data  in  file  goes  from  north  to  south 
zrev  -  z  data  in  file  goes  from  top  to  bottom 
sequential  -  FORTRAN  unformatted  sequential 
template  -  file  name  templates  in  use 
big_endian  -  data  is  big_endian 
little_endian  -  data  is  little_endian 

Sequential:  The  sequential  option  may  be  set  if  the  data  were 
written  using  unformatted  (f77)  I/O. 

Template  (for  multiple  file  time  series):  GrADS  allows  you  to 
handle  many  actual  data  files  as  one  GrADS  file,  if  the  individual 
data  files  are  in  a  GrADS  readable  format,  and  if  the  files  are 
split  along  time.  In  the  initial  implementation,  the  time(s)  that 
are  in  each  file  are  indicated  by  the  file  name. 

An  example  of  this  might  be  hourly  data,  where  each  24  hours 
has  been  placed  in  a  separate  file.  Each  file  is  named 
this  way: 

lmay92.dat 

2may92.dat 

etc. 

You  indicate  to  GrADS  that  there  are  multiple  files  in  this 
time  series  by  giving  a  substitution  template  as  the  file  name: 

dset  %dl%mc%y2.dat 

and  giving  an  options  record  that  looks  like: 
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options  template 

and  specifying  the  time  range  and  increment  in  the 
tdef  record; 

tdef  72  linear  0zlmayl993  Ihr 

GrADS  will  figure  out  automatically  that  there  are  24  times 
in  each  file,  and  what  file  names  correspond  to  what  times. 

As  you  display  data,  GrADS  will  only  open  one  file  at  a  time. 

As  you  change  times  such  that  another  file  is  referred  to, 
the  open  file  is  closed,  and  the  new  file  is  opened. 

Valid  substitutions  are: 

%y2  -  2  digit  year  (last  2  digits) 

%y4  -  4  digit  year 

%ml  -  1  or  2  digit  month 

%m2  -  2  digit  month  (leading  zero  if  needed) 

%mc  -  3  character  month  abbreviation 

%dl  -  lor  2  digit  day 

%d2  -  2  digit  day 

%hl  -  1  or  2  digit  hour 

%h2  -  2  digit  hour 

This  support  works  on  all  supported  GrADS  data  types  (GrADS 
gridded,  GRIB,  GrADS  station  data).  If  you  specify  file  format 
options,  the  options  must  apply  equally  to  each  file. 

Byte  ordering:  You  may  specify  the  actual  byte  ordering  of  the  binary  data  file: 

options  big_endian 
or 

options  little_endian 

Indicates  the  binary  data  file  is  in  reverse  byte  order  from  the  normal  byte  order 
of  the  machine. 

If  the  data  is  already  in  the  correct  order,  no  conversion 
is  performed.  If  you  then  move  the  data  set  to  another  machine, 
conversion  will  be  automatically  performed.  This  would  happen 
if  you  sent  a  file  in  binary  format  from,  for  example,  a  Sun  to  a  PC. 

These  options  were  added  to  facilitate  moving  data  files  and  descriptor 
files  between  machines.  Putting  this  keyword  in  the  descriptor  file  tells 
GrADS  to  correct  the  byte  order  as  the  data  is  being  read. 
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XDEF  number  <LINEAR  start  increment> 

<LEVELS  value-list> 

Defines  the  mapping  between  grid  values  and  longitude. 

Specifically; 

number  -  the  number  of  grid  values  in  the  X  direction, 
specified  as  an  integer  number. 

Must  be  >=  1. 

LINEAR  or  LEVELS  —  Indicates  the  grid  mapping  type. 

For  LINEAR: 

start  —  the  starting  longitude,  or  the  longitude  for  X  =  1. 
Specified  as  a  floating  point  value,  where  negative 
indicates  degrees  west. 

increment  --  the  spacing  between  grid  value  in  the  X  direction. 

It  is  assumed  that  the  X  dimension  values  go  from 
west  to  east.  Specified  as  a  positive  floating  value. 

For  LEVELS: 

value-list  —  List  of  'number'  values  representing  the 

longitude  of  each  X  dimension.  May  start  and 
continue  on  the  next  record  in  the  descriptor 
file  (records  may  not  be  >  80  characters). 

There  must  be  at  least  2  levels  (otherwise 
use  LINEAR  mapping). 

YDEF  number  mapping  start  <increment> 

<LEVELS  value-list> 

Defines  the  mapping  between  grid  values  and  latitude. 

Specifically: 

number  —  the  number  of  grid  values  in  the  X  direction, 
specified  as  an  integer  number. 

mapping  —  mapping  type,  specified  as  a  keyword. 

Valid  are: 

LINEAR  -  Linear  mapping 
GAUSR15  —  Gaussian  R15  latitudes 

GAUSR20  —  Gaussian  R20  latitudes 

GAUSR30  —  Gaussian  R30  latitudes 

GAUSR40  —  Gaussian  R40  latitudes 

start  —  For  LINEAR  mapping,  the  starting  latitude,  ie  the 
latitude  for  Y  =  1,  and  is  specified  as  a  floating 
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point  value,  with  negative  indicating  degrees  south. 

For  GAUSRxx  mapping,  the  start  value  indicates 
the  first  gaussian  grid  number,  where  1  would  be  the 
southernmost  gaussian  grid  latitude, 
increment  ~  the  spacing  between  grid  values  in  the  Y  direction. 
It  is  assumed  that  the  Y  dimension  values  go  from 
south  to  north.  Specified  as  a  positive  floating 
point  value.  Used  only  for  LINEAR  mapping. 

For  LEVELS: 

value-list  -  List  of  'number'  values  representing  the 
latitude  of  each  X  dimension.  May  start  and 
continue  on  the  next  record  in  the  descriptor 
file  (records  may  not  be  >  80  characters). 

There  must  be  at  least  2  levels  (otherwise 
use  LINEAR  mapping). 

Examples  of  specifying  GAUSRxx  mapping: 

YDEF  20  GAUSR40  15 

Indicates  that  there  are  20  Y  dimension  values  which  start 
at  Gaussian  Latitude  15  (64.10  south)  on  the  Gaussian  R40  grid. 
Thus  the  20  values  would  correspond  to  Latitudes: 

-64.10,  -62.34,  -60.58,  -58.83,  -57.07,  -55.32,  -53.56, 

-51.80,  -50.05,  -48.29,  -46.54,  -44.78,  -43.02,  -41.27, 

-39.51,  -37.76,  -36.00,  -34.24,  -32.49,  -30.73 

YDEF  102  GAUSR40  1 

The  entire  gaussian  grid  is  present,  starting  at  the  southernmost 
latitude  (-88.66). 

ZDEF  number  mapping  <start  increment> 

<value-list> 

Defines  the  mapping  between  grid  values  and  pressure  level. 
Specifically: 

number  -  the  number  of  grid  values  in  the  X  direction, 
specified  as  an  integer  number, 
mapping  -  mapping  type,  specified  as  a  keyword. 

Valid  are: 

LINEAR  —  Linear  mapping 

LEVELS  —  Arbitrary  pressure  levels 
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start  -  when  mapping  is  LINEAR,  this  is  the  starting  value,  or 
the  value  when  Z=l. 

increment  —  when  mapping  is  LINEAR,  the  increment  in  the  Z 
direction,  or  from  lower  to  higher.  This  may  be  a 
negative  value,  for  example: 

ZDEF  10  LINEAR  1000  -100 

indicating  that  the  data  is  for  levels  1000,  900,  800, 

700,  etc. 

value-list  —  when  the  mapping  is  LEVELS,  the  specific  levels 
are  simply  listed  in  ascending  order.  If  there 
is  only  one  level,  use  LINEAR,  since  LEVELS 
implies  at  least  two  levels. 

TDEF  number  LINEAR  start-time  increment 

Defines  the  mapping  between  grid  values  and  time. 

Specifically: 

number  -  the  number  of  times  in  the  data  set. 

Specified  as  an  integer  number, 
start-time  -  The  starting  date/time  value,  specified  in  GrADS 
absolute  date/time  format.  This  is  the  value 
when  T=l.  The  date/time  format  is: 

hh:mmZddmmmyyyy 

where: 

hh  =  hour  (two  digit  integer) 
mm  =  minutes  (two  digit  integer) 
dd  =  day  (one  or  two  digit  integer) 
mmm  =  month  (jan,  feb,  mar,  apr,  may,  jun,  jul, 
aug,  sep,  oct,  nov,  dec) 
yyyy  =  year  (two  or  four  digit  integer. 

two  digits  implies  a  year  between 
1950  and  2049). 

If  not  specified,  hh  defaults  to  00,  mm  defaults  to  00, 
and  dd  defaults  to  1.  The  month  and  year  must  be 
specified.  No  intervening  blanks  are  allowed  in  a 
GrADS  absolute  date/time. 

Examples: 
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12Z1JAN1990 

14:20Z22JAN1987 

JUN1960 

increment  —  time  increment.  Specified  in  GrADS  time  increment 
format: 

vvkk 

where: 

vv  =  an  integer  number,  1  or  2  digits 
kk  =  an  increment  keyword, 
mn  =  minutes 
hr  =  hours 
dy  =  days 
mo  =  months 
yr  =  year 

Examples: 

20mn  —  increment  is  20  minutes 
Imo  —  increment  is  1  month 
2dy  —  increment  is  2  days 

Some  examples  of  a  TDEF  statement: 

TDEF  24  LINEAR  00Z01JUN1987  IHR 

The  data  set  has  24  times,  starting  at  OOZ  on  1  Jun,  1987, 
with  an  increment  of  1  hour. 

TDEF  30  LINEAR  2JUN1988  IDY 

The  data  set  has  30  times,  starting  at  OOZ  on  2  Jun,  1988, 
with  an  increment  of  1  day. 
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VARS  number 


Indicates  the  start  of  the  records  describing  the  variables  in  the  data  set. 
number  —  the  number  of  variable  records 
Each  variable  record  is  in  the  following  format; 
abrev  levs  units  description 

abrev  —  a  1  to  8  character  abbreviation  for  this  variable. 

This  abbreviation  must  start  with  an  alphabetic 
character  (a-z)  and  be  composed  of  alphabetic 
characters  and  numbers.  This  abbreviation  will  be 
the  "name"  the  variable  is  accessed  by  from  within 
GrADS. 

levs  —  an  integer  value  specifying  the  number  of  levels 
this  variable  has  in  the  data  set.  It  may  not 
exceed  the  number  of  levels  in  the  ZDEF  statement. 

A  levs  value  of  0  indicates  this  variable  has  one 
"level"  that  does  not  correspond  to  a  vertical  level. 

An  example  would  be  a  surface  variable. 

units  —  Reserved  for  future  use.  Put  a  value  of  99  here. 

description  -  A  text  description  of  the  variable,  max  40  characters. 

After  the  last  variable  record  comes  the  ENDVARS  statement.  This  ends  the  GrADS 
data  descriptor  file. 

GRIB  data  support:  GrADS  will  now  read  some  GRIB  data  files  directly, 
particularly  those  being  produced  by  NMC.  The  GRIB  records  must 
contain  lat-lon  grids.  A  detailed  discussion  on  how  to  use  GRIB  data  files  is 
provided  in  GrADS  documentation  titled:  update.  151. 

A.2.  Station  Data  Sets. 

Station  data  sets  are  written  to  a  binary  file  one  report  at  time.  The  only  ordering 
required  is  that  the  station  reports  be  grouped  within  the  file  into  some  time  interval. 
For  example,  the  time  interval  for  upper  air  observations  might  be  12  hours. 

Variables  within  each  report  are  spht  into  two  groupings.  Each  variable  is  either  a 
surface  variable,  thus  can  be  reported  at  most  once  per  report,  or  it  is  a  level  dependent 
variable,  thus  can  be  reported  at  a  number  of  different  levels  within  one  report. 

The  format  of  a  station  report  in  the  binary  station  data  file  is: 
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-  A  header  which  provides  information  about  the 
location  of  the  station. 

-  Surface  variables,  if  any 

-  Level  dependent  variables,  if  any 


The  header  is  described  by  the  following  C  language  data  structure; 


struct  rpthdr  { 
char  id[8]; 
float  lat; 
float  Ion; 
float  t; 
int  nlev; 
int  flag; 

}; 


/*  Character  station  id  */ 

/*  Latitude  of  report  */ 

/*  Longitude  of  report  */ 

/*  Time  in  relative  grid  units  */ 

/*  Number  of  levels  following  */ 
/*  Level  independent  var  set  flag  */ 


A  detailed  description  of  each  header  entry  follows: 

id  -  The  station  identifier.  This  is  a  1  to  7  character 
identifier  that  should  identify  the  station  uniquely. 

It  may  be  assigned  arbitrarily;  ie.  the  stations  could 
be  numbered  in  some  arbitrary  order, 
lat  -  The  Y  dimension  location  of  the  station  in  world 
coordinates,  typically  latitude. 

Ion  -  The  X  dimension  location  of  the  station  in  world 
coordinates,  typically  longitude, 
t  -  The  time  of  this  report,  in  relative  grid  units. 

This  refers  to  the  way  the  stations  are  grouped  in 
time.  For  example,  if  you  are  working  with  surface 
airways  reports,  you  would  probably  have  a  time  grouping 
interval  of  one  hour.  If  you  wanted  to  treat  the 
report  times  of  each  report  as  being  exactly  on  the 
hour,  you  would  set  t  to  0.0.  If  the  report  was 
for  12:15pm,  and  you  were  writing  the  time  group  for 
12pm,  you  would  set  t  to  be  0.25.  Thus,  t  would 
typically  have  the  range  of  -0.5  to  0.5. 
nlev  -  Number  of  data  groups  following  the  header.  This 
is  the  count  of  the  one  surface  group,  if  present, 
plus  the  number  of  level  dependent  groups. 

Is  set  to  zero  to  mark  the  end  of  a  time  group  in  the 
file. 

flag  -  If  zero,  there  are  no  surface  variables  following  the 
header.  If  one,  then  there  are  surface  variables 
following  the  header. 
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Following  the  header,  the  data  for  this  report  is  written.  The  first  group  of  data  would 
be  all  the  surface  variables  if  present.  Whether  or  not  the  surface  variable  (if  any)  are 
present  is  determined  by  the  flag  in  the  header.  If  present,  then  all  the  surface 
variables  must  be  written  -  missing  variables  should  have  the  missing  data  value 
provided.  Thus,  each  surface  variable  group  will  be  the  same  size  for  each  report  in 
the  file. 

The  surface  variables  are  written  out  as  floating  point  numbers.  The  ordering  of  the 
variables  must  be  the  same  in  each  report,  and  is  the  ordering  that  will  be  given  in  the 
data  descriptor  file. 

Following  the  surface  variable  group,  any  number  of  level  dependent  groups  may  be 
written.  The  number  of  total  data  groups  is  provided  in  the  header.  Each  level 
dependent  group  must  have  all  the  level  dependent  variables  present,  even  if  they  are 
filled  with  the  missing  data  value.  Thus,  each  level  dependent  group  will  be  the  same 
size  for  all  levels  and  all  reports  in  the  file. 

The  level  dependent  group  is  written  out  as  follows: 

level  —  floating  point  value  giving  the  Z  dimension 
value  in  world  coordinates  for  this  level, 
variables  —  The  level  dependent  variables  for  this  level. 

After  all  the  reports  for  one  time  grouping  have  been  written,  a  special  header  (with  no 
data  groups)  is  written  to  indicate  the  end  of  the  time  group.  The  header  has  an  nlev 
value  of  zero.  The  next  time  group  may  then  start  immediately  after.  A  time  group 
with  no  reports  would  still  contain  the  time  group  terminator  header  record  (ie,  two 
terminators  in  a  row). 

GrADS  station  data  files  must  be  written  as  UNIX  stream  data  sets  without  any 
imbedded  record  descriptor  information.  This  is  easily  done  from  a  C  program.  From 
a  FORTRAN  program,  it  usually  requires  a  system-dependent  option  in  the  OPEN 
statement.  For  example,  in  DEC  FORTRAN  one  can  use  the 

RECORDTYPE=’STREAM' 

option  to  avoid  having  record  descriptor  information  imbedded  in  the  output  file. 
Examples  of  C  and  FORTRAN  programs  to  create  station  data  sets  are  provided  later 
in  this  document. 


Station  Data  Descriptor  File 

The  format  for  the  data  descriptor  file  for  station  data  is  similar  to  the  format  for  a 
gridded  data  set.  An  example  of  a  station  data  descriptor  file  is: 


dset  ^ua.reps 
dtype  station 
stnmap  ^ua.map 
undef  -999.0 

title  Real  Time  Upper  air  obs 

tdef  10  linear  12zl8janl992  12hr 

vars  12 

sip  0  99  SLP 

ts  0  99  Temps 

ds  0  99  Dewpoints 

us  0  99  U  Winds 

vs  0  99  V  Winds 

z  1  99  Heights 

t  1  99  Temps 

d  1  99  Dewpoints 

u  1  99  U  Winds 

V  1  99  V  Winds 

endvars 

Note  the  differences  between  this  descriptor  file  and  a  grid  descriptor  file: 

DTYPE  record  -  specify  a  data  type  of:  station. 

STNMAP  record  —  gives  the  file  name  of  the  station  mapping 
file.  This  file  is  created  by  the 
stnmap  utility,  which  will  be  described 
later. 

XDEF,  YDEF,  ZDEF  records  -  not  specified. 

TDEF  record  —  describes  the  time  grouping  interval  and  the 
number  of  time  groups  in  the  file. 

VAR  records  —  surface  variables  must  come  first,  and 
are  given  a  zero  for  the  number-of-levels 
field.  Level  dependent  variables  are  hsted 
after  the  surface  variables,  and  are  given 
a  one  in  the  number-of-levels  field. 

Byte-ordering  control  for  station  data  files:  You  may  specify 
byte  ordering  (big_endian,  or  little_endian)  for  station  data  files 
using  the  OPTIONS  record.  The  stnmap  utility,  and  GrADS, 
will  perform  the  necessary  conversion.  Station  map  files  must 
be  created  on  the  machine  where  they  are  to  be  used. 

STNMAP  Utility 

Once  the  data  set  has  been  written,  and  the  descriptor  file  created,  you  should  then 
create  the  station  map  file  by  running  the  stnmap  utility.  Yhis  utility  writes  out  hash 
table  and/or  link  hst  information  that  allows  GrADS  to  access  the  report  data  more 
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efficiently.  The  utility  will  prompt  for  the  name  of  the  data  descriptor  file. 

If  you  change  the  data  file  —  perhaps  by  appending  another  time  group  —  you  will 
also  have  to  change  the  descriptor  file  to  reflect  the  changes  -  the  new  number  of 
times  for  example  -  and  then  rerun  the  stnmap  utility. 


A.3.  Examples  of  Creating  a  Gridded  Data  Set. 

On  a  workstation,  the  binary  GrADS  data  sets  need  to  be  created  as  a  'stream'  data  set, 
ie,  it  should  not  have  the  normal  FORTRAN  record  descriptor  words  imbedded  in  it. 
This  can  be  done  from  FORTRAN  using  direct  access  FO: 

REAL  Z(72,46,16) 


OPEN(8,nLE='grads.dat',FORM='UNFORMATTED', 
&  ACCESS='DIRECT',RECL=72*46) 


IREC=1 
DO  101=1,18 

WRITE  (8,REC=IREC)  ((Z(J,K,I),J=1,72),K=1,46) 

IREC=IREC+1 
10  CONTINUE 

This  example  writes  out  16  levels  of  one  variable  to  a  file  in  direct  access  format.  We 
are  really  writing  the  data  out  sequentially,  and  using  direct  access  to  avoid  having  the 
record  descriptor  words  written.  There  may  be  options  in  your  compiler  to  do  this 
more  directly,  or  you  may  wish  to  write  the  data  using  a  C  program. 


Another  simple  sample  might  be: 

REALX(IOO) 

DO  10 1=1,100 
X(I)=I 

10  CONTINUE 

OPEN(8,FILE='samp.dat',FORM='UNFORMATTED',ACCESS='DIRECT', 
&  RECL=100) 

WRITE  (8,REC=1)  X 

STOP 

END 


The  associated  descriptor  file: 
DSET  samp.dat 
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TITLE  Sample  Data  Set 

UNDEF  -9.99E33 

XDEF  100  LINEAR  1  1 

YDEF  1  LINEAR  1  1 

ZDEF  1  LINEAR  1  1 

TDEF  1  LINEAR  1JAN2000  IDY 

VARS  1 

X  0  99  100  Data  Points 
ENDVARS 

Once  created,  you  can  use  this  data  set  to  experiment  with 
GrADS  data  functions,  such  as: 

display  sin(x/50) 


A.4.  Examples  of  Creating  Station  Data  Sets. 

Lets  say  you  have  a  data  set  with  monthly  rainfall: 


Year  Month  Stid 

Lat  Lon 

Rainfall 

1980 

1  QQQ 

34.3  -85.5 

123.3 

1980 

1  RRR 

44.2  -84.5 

87.1 

1980 

1  SSS 

22.4  -83.5 

412.8 

1980 

1  TTT 

33.4  -82.5 

23.3 

1980 

2  QQQ 

34.3  -85.5 

145.1 

1980 

2  RRR 

44.2  -84.5 

871.4 

1980 

2  SSS 

22.4  -83.5 

223.1 

1980 

2  TTT 

33.4  -82.5 

45.5 

A  FORTRAN  program  in  DEC  FORTRAN  to  write  this  data  set  in  GrADS  format 
might  be: 


CHARACTER'S  STID 
C 

OPEN  (8,NAME='rain.ch’) 

OPEN(10,NAME='rain.dat',FORM='UNFORMATrED', 
&  RECORDTYPE='STREAM') 


IFLAG  =  0 


c 

C  Read  and  Write 
C 

10  READ  (8,9000,END=90)  IYEAR,IMONTH,STID,RLAT,RLON,RVAL 
9000  FORMAT  (I4,3X,I2,2X,A8,3F8.1) 

IF  (IFLAGEQ.O)  THEN 
IFLAG  =  1 
lYROLD  =  lYEAR 
IMNOLD  =  IMONTH 
ENDIF 
C 

C  If  new  time  group,  write  time  group  terminator. 

C  Assuming  no  empty  time  groups. 

C 

IF  (IYROLD.NE.IYEAR.OR.IMNOLD.NE.IMONTH)  THEN 
NLEV  =  0 

WRITE  (10)  STID,RLAT,RLON,TIM,NLEV,NFLAG 

ENDIF 

lYROLD  =  lYEAR 
IMNOLD  =  IMONTH 
C 

C  Write  this  report 
C 

TIM  =  0.0 
NLEV  =  1 
NFLAG  =  1 

WRITE  (10)  STID,RLAT,RLON,TIM,NLEV,NFLAG 
WRITE  (10)  RVAL 
GO  TO  10 
C 

C  On  end  of  file  write  last  time  group  terminator. 

C 

90  CONTINUE 
NLEV  =  0 

WRITE  (10)  STID,RLAT,RLON,T[M,NLEV,NFLAG 

STOP 

END 

For  a  different  compiler,  you  would  need  to  determine  the  appropriate  OPEN 
statement  to  write  a  stream  data  set. 


An  equivalent  C  program  might  be: 
#include  <stdio.h> 
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/*  Structure  that  describes  a  report  header  in  a  stn  file  */ 


struct  rpthdr  { 
char  id[8]; 
float  lat; 
float  Ion; 
float  t; 
int  nlev; 
int  flag; 

}  hdr; 


/*  Character  station  id  *! 

/*  Latitude  of  report  */ 

/*  Longitude  of  report  */ 

/*  Time  in  relative  grid  units  */ 

/*  Number  of  levels  following  */ 
/*  Level  independent  var  set  flag  */ 


main  ()  { 

FILE  *ifile,  *ofile; 
char  rec[80]; 

int  flag,year, month, yrsav,mnsav,i; 
float  val; 


/*  Open  files  */ 

ifile  =  fopen  ("rain.ch","r"); 
ofile  =  fopen  ("rain.dat", "wb"); 
if  (ifile==NULL  II  ofile==NULL)  { 
printf  ("Error  opening  files\n"); 
return; 

} 

/*  Read,  write  loop  */ 


flag  =  1; 

while  (fgets(rec,79,ifile)!=NULL)  { 
/*  Format  conversion  */ 


sscanf  (rec,"%i  %i  ",&year,&month); 

sscanf  (rec+20,"  %g  %g  %g",&hdr.lat,&hdr.lon,&val); 

for  (i=0;  i<8;  i++)  hdr.id[i]  =  rec[i+ll]; 


/*  Time  group  terminator  if  needed  */ 


if  (flag)  { 
yrsav  =  year; 
mnsav  =  month; 
flag  =  0; 

} 


if  (yrsav  !=year  II  mnsav  !=month) 
hdr.nlev  =  0; 


{ 
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fwrite  (&hdr,sizeof(stmct  rpthdr),  1,  ofile); 

} 

yrsav  =  year; 
mnsav  =  month; 

/*  Write  this  report  */ 

hdr.nlev  =  1 ; 
hdr.flag  =  1 ; 
hdr.t  =  0.0; 

fwrite  (&hdr,sizeof(struct  rpthdr),  1,  ofile); 
fwrite  (&val,sizeof(float),  1,  ofile); 

} 

hdr.nlev  =  0; 

fwrite  (&hdr,sizeof(stract  rpthdr),  1,  ofile); 

} 

Once  the  binary  data  file  has  been  written,  create  the  descriptor  file.  It  would  look 
something  Like  this: 

dset  rain.dat 
dtype  station 
stnmap  rain.map 
undef  -999.0 
title  Rainfall 

tdef  12  linear  janl9801mo 
vars  1 

p  0  99  Rainfall 
endvars 

Then  run  the  stnmap  utility  to  create  the  station  map  file.  You  can  then  open  and 
display  this  data  from  within  GrADS. 

A.5.  Examples  of  ‘stations.ua’  and  ‘stations.sfc’  Data  Sets. 

The  upper  air  station  file  is  an  ASCII  list  that  includes  the  Stid,  Lat,  Lon,  and  station 
title  as  shown  here: 

NZY  32.7  -1 17.2  North  Island 
MIR  32.9-117.1  San  Diego 


EDW  34.9 -117.9  Edwards  AFB 
LAS  36.6 -116.0  Las  Vegas 
OAK  37.7 -122.2  Oakland 
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The  surface  station  file  is  also  an  ASCII  list,  but  only  includes  the  Stid,  Lat,  and  Lon  as 
shown  here: 


ACV  41.0-124.1 
BAB  39.1  -121.4 


UKI  39.1  -123.2 
WMC  40.9-117.8 
YUM  32.7-114.6 


A-21 


This  Page  Intentionally  Left  Blank. 
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Appendix  B:  NEBS  Scripts 


NEBS.GS  . B-3 

OUTPUT.GS . B-35 


This  Page  Intentionally  Left  Blank. 
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Ions  =  Jons  %  ' '  %  siibwrd(result,4) 
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